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SECTION  I 


INTRODUCTION 


For  enhancing  verticahy-incident  P-waves,  the  application 
of  multichannel  filter  techniques  to  seismometer  systems  consisting  of 
both  horizontal-  and  vertical-component  instruments  presents  an  attrac¬ 
tive  alternative  to  the  use  of  systems  consisting  solely  of  vertical  seis¬ 
mometers.  The  advantage  of  using  horizontals  is  that  it  permits  s  type  of 
noise-reduction  processing  which  (1)  exploits  coupling  between  horizontal 
and  vertical  components  of  ambient  seismic  noise  and  (2)  preserves  the 
exact  form  of  the  P-wave  signal. 

In  this  report  we  shall  discuss  the  theoretical  capabilities 

of  systems  consisting  of  a  single  vertical  seismometer  v  and  one  or  more 

horizontal  seismometers  h,,  h_,  .  .  .  h_  We  shall  restrict  our  atten- 

12  M 

tion  to  the  following  type  of  processings  we  apply  to  each  horizontal  seis¬ 
mometer  output  h.  a  linear  time-invariant  filter  g.,  and  subtract  all  of  the 

J  J 

filter  outputs  g.  ®  h.  from  v  to  obtain  the  output  o  of  our  processor.  Thus, 

o  =  v  -  g1  ®  hj  -  g2  ®  h2  -  .  .  .  gM  ®  hM  (Figure  1.  1). 

Evidently,  this  processing  scheme  leads  to  no  degradation  of  the  P-wave 
signal  since  the  vertically-incident  P-wave  has  no  horizontal  component. 

In  order  to  achieve  maximum  noise  suppression,  the  filters  g^  must  be 
designed  in  such  a  way  that  (in  the  absence  of  signal) 

T 

/  H‘>f  * 

-.T 

is  minimized.  Such  optimal  filters  are  called  Wiener  prediction  filters; 
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they  are  determined  by  the  matrix  equation  (Burg,  ’964) 
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where  the  Gj  are  the  frequency-domain  responses  of  the  filters  g.,  and 
the  $  are  the  auto-  and  crosspower  spectra  of  the  sei.  mometer 

traces. 


The  simplest  processor  of  this  type  is  the  2 -component 
seismometer  at  point  location;  i.  e.  ,  a  single  vertical  seismometer  v  and 
a  single  horizontal  seismometer  h  located  at  the  same  point  on  the  earth's 
surface.  Here,  the  output  of  the  processor  is  simply 


o  =  v  -  g  ©  h 
where  g(t)  is  determined  by 


G(l)  =  »yh/ » 


hh 


Part  I  of  this  report  is  devoted  to  a  detailed  consideration  of  the  capabili¬ 
ties  and  limitations  of  this  2 -channel  processor  in  ambient  trapped-mode 
noise.  The  feasibility  of  using  a  2 -component  seismometer  as  a  P-wave 
enhancer  is  suggested  by  the  observation  that  complete  suppression  of 
noise  would  be  achieved  with  this  system  in  the  ideal  situation  thav  the 
noise  is  simply  a  unidirectional  single-mode  Rayleigh  wave.  In  that 
case,  v  would  be  purely  a  convolution  of  h;  that  is,  v  -  k  ©  h  for  some 
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function  k(l),  and  taking  g  -  k,  we  get  o(t)  =  U.  One  should  expect  that 
deviation  from  this  ideal  performance  would  be  due  principally  to  the 
following  factors: 

(1)  If  the  Rayleigh  wave  field  is  not  unidirectional,  but  comes  from 
a  distribution  of  azimuths,  there  will  be  less  than  perfect  perform¬ 
ance  even  if  all  the  noise  propagates  in  a  single  mode.  In  fact,  in  iso¬ 
tropic  noise,  uncorrelated  by  azimuth,  no  noise  reduction  whatsoever  can 
be  achieved  with  the  system. 

(2)  Realistically,  v  and  h  must  be  assumed  to  contain  components 
of  random  noise  uncorrelated  with  the  Rayleigh  waves  (e.  g.  ,  amplifier 
noise). 

(2)  The  presence  on  the  horizontal  of  Love  (SH)  waves,  statistically 
independent  from  the  Rayleigh  waves,  reduces  performance.  (One  may 
conceive  of  Love  waves  on  h  as  a  special  case  of  (2),  but  it  will  be  seen 
that  it  may  sometimes  be  misleading  to  do  so.) 

(4)  The  amplitude  gains  of  seismometers  vary  slowly  with  time,*  this 
places  a  limitation  on  the  effectiveness  of  a  system  in  which  the  prediction 
filter,  once  it  has  been  designed,  is  left  online  for  an  extended  period  of 
time. 

(5)  Several  modes  of  Rayleigh-type  noise  may  be  present,  each  having 
a  different  transfer-function  relating  the  vertical  and  horizontal  component. 

The  quantitative  effects  of  (1)  through  (5)  on  the  prediction 
capability  of  the  2 -component  seismometer  at  point  location  are  discussed 
in  Sections  III,  IV,  V,  VI,  and  VII,  respectively.  Section  II  is  devoted  to 
a  summary  of  communication  theoiy  formulae  and  notation  used  in  the  rest 
oi  the  paper. 
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Part  II  begins  in  Section  IX  with  a  derivation  of  general 

formulae  for  the  response,  in  single-mode  noise,  of  a  planar  array  of 

seismometers  v,  h  ,  h. ,  .  .  .,  h.  .,  where  the  horizontals  h.  may  have 

1  2  M  j 

arbitrary  orientations  and  may  be  separated  from  each  other  or  from 
the  vertical  v.  Certain  arrays  of  this  type  will  be  shown  to  give  good 
Performance  even  in  isotropic  noise.  In  Section  X,  the  formulae  from 
Section  IX  are  used  to  calculate  the  response  in  isotropic  single-mode 
noise  (with  the  possible  addition  of  Love  waves)  of  arrays  consisting  of 
radially-oriented  horizontal  seismometers  regularly  spaced  around  a 
ring  having  the  vertical  seismometer  at  its  center.  It  will  be  shown  that 
such  array  geometries  are  optimal  in  a  certain  sense. 

The  sections  of  this  report  are,  for  the  most  part,  inde¬ 
pendent  of  each  other  with  the  following  exceptions:  Section  II  is  intended 
for  reference  only;  Section  III  should  be  read  first;  Section  IV  should  be 
read  before  Section  V;  and  Section  IX  should  be  read  before  Section  X. 
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SECTION  II 


COMMUNICATION  THEORY  FORMULAS  AND  NOTATION 

Primarily  for  the  purpose  of  fixing  notation,  we  shall  state 
several  formulas  to  be  used  in  the  rest  of  this  paper. 

A.  CORRELATION  FUNCTIONS 

Let  x(t)  be  a  stationary  random  process.  Then  the  auto¬ 
correlation  function  rc  (t)  of  x(t)  is  defined  as 

xx 


T 

xx(t)  -  T™  /  *(T,x(T-t)dT 


By  ergodicity,  the  autocorrelation  function  rc  (t)  depends  only  upon  the 


xx 


ensemble  to  which  x(t)  belongs. 


The  autopower  spectrum  $^(f)  of  x(t)  is  the  Fourier  trans 


form  of  cp  (t); 
xx 


*  -<f)  =  f  ?  (t)e'i2nf‘dt 

k.  J  XX 

—  00 


XX 


Similarly,  if  x(t)  and  y(t)  are  stationary  random  processes, 

the  crosscorrelation  function  ro  (t)  of  x  and  y  is 
-  xy 


(t)  -  ^  r x(T)y  (T_t) 

-T 


xy 


dT 


and  the  crosspower  spectrum  $  ^(f)  nf  x  and  y  is 


$  (f)  =  f 
xy 


„A  -iZrrft  . 

T  (t)e  dt 

xy 


(2.  0 


(2.Z) 


(2.  3) 


(2.4) 
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The  crosspower  spectrum  is  easily  seen  to  satisfy 


and 


Vi+V2’ 


(f)  =  a, 


$  (f) 

xly 


+  $ 


V 


»  (£)  =  ♦  (£) 

yx  xy 


denotes  complex 
conjugation 


B.  LINEAR  FILTERS 

The  output  z(t)  of  the  most,  general  linear  filter  applied  to 
x(t)  may  be  expressed  as 

r 00 

z(t)  =  (k®x)(t)  =  J  k(T)x(t  -^)  dT 

•  CO 

where  k(T)  is  a  function  such  that 

I  |k(T )  |  dT  <  ® 

J  m 


The  Fourier  transform  K(f)  of  k(T)  is  called  the  transfer 
function  of  the  filter  k. 


C.  CROSSPOWER  SPECTRA  OF  FILTERED  TIME  FUNCTIONS 


Suppose  x(t)  and  y(t)  are  stationary  random  processes,  and 
g(t)  and  k(t)  are  linear  filters.  Let  G(f)  and  K(f)  be  the  transfer  functions 
corresponding  to  g  and  k  respectively.  Then  it  is  easily  verified  from 
(2.  1)  -  (2.  6)  that 


and 


$x,k© 
©  x, 

©  x, 


y(f)  =  K*  (f)  *xy(f) 
y(f)  =  G(f)  $xy(f) 

k%(£)  =  G(f)  K"  (f)  $xy(f) 
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(2.5) 


(2.  6) 


(2.7) 


In  particular, 


g  v35  x,  g 


<f)  =  1  G(f ) 


xx 


M) 


(2.8) 


D.  COHERENCE 


coherence 


C 

xy 


Let  x  (t)  and  y(t)  be  stationary  random  process 
(f)  between  x  and  y  is  defined  as  follows: 


es. 


The 


C 

xy 


(f) 
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-  JS1 
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xx 


(f)  *  (f) 
_ Z2£_ 

(f)  *  (f) 

yy 
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$  (f)  $ 

xx  yy 


(f) 


(2.9) 


It  may  be  shown  that  C(f)  satisfies  0  £  C(f)  £  1  for  all  real  f. 

If  g  and  k  are  linear  filters,  it  is  immediate  from  (2.  7)  and 

(2.  9)  that 


g©x,  k©y(f) 


C  (f) 
xy 


(2.  10) 


That  is,  the  coherence  between  x  and  y  is  invariant  with  respect  to 
transformation  of  x  and  y  by  arbitrary  linear  filters. 

E.  PREDICTION 

Let  v(t)  and  h(t)  be  stationary  random  processes  (e.g.,  vertical 
and  horizontal  components  of  ambient  seismic  noise  at  some  location).  We 
wish  to  find  the  best  frequency  filter  g(t)  to  apply  to  h(t)  such  that  the  filtered 
signal  li(t)  =  (g©h)(t)  minimizes  the  expression 


lim 

■p-4  00 
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If  we  let  e(t)  =  v(c)  -  n(t),  and  G(f)  =  Fourier  transform  of  g(t),  our  problem 
is  equivalent  to  finding  G(f)  such  that 


„  Qp 


is  minimized. 


(2.11) 


Using  (2.  5),  we  have 


$  (f)  =  t  ,/  ,»(f) 

ee  v  -h  ,  v-h 


(2.  12) 


-  *vv<£>  -  *vh'<£>  -  Vv<£>  +  Vh'(£> 

-  -  2Re  (#h'v<£>)  +  Vh'  » 


Hence  by  (2.  7) , 


W£)  *  #VV(£)  -  2Re  (G(f)  *hv<1))  +  lG(f)!2\h(f) 


(2.13) 


The  last  expression,  the  integrand  in  (2.  11),  is  non-negative 
for  all  f.  Hence  to  minimize  (2.  11),  we  minimize  (2.  13)  at  each  f.  Letting 

G(f)  =  cd  +  i  lit 


in  (2.  13)  and  setting 


'(*ee<£>)  »(*..«) 


=  0 


we  find 


cp  = 


(*hv<£>) 


_  Im 


W£> 


.  t  « 


(*hv<£>) 


*hh<£> 


Hence  the  optimum  prediction  filter  is  given  by 


G(f)  = 


*vh<f) 

Vf> 


(2.  14) 


Substitution  of  (2.  14)  into  (2.  13)  shows  that  the  autopower 
spectrum  of  the  least-mean-square  error  is 


minimal  $  (f) 

ee 


*  (f) 

$  (f)  i .  — - — — 

VV  $vv(f)  $hh(f) 


(2.15) 


=  $  (f) 

VV 


|^1  -  C(f)J 


where  C(f)  =  coherence  between  h  and  v. 


Because  of  (2.  15),  we  call  I  -C(f)  the  relative  prediction 
error  for  the  frequency  f. 
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SECTION  III 


a 

SINGLE  MODE  RAYLEIGH  WAVE  FIELD"' 


Let  O  be  the  origin  of  an  X-Y  coordinate  system  in  the  plane 
of  the  earth's  surface  and  let  Z  be  a  vertical  axis  through  O,  positive  down¬ 
ward.  If  Q  is  any  point  in  the  X-Y  plane,  the  azimuth  of  Q  is  the  angle  0 
measured  counterclockwise  from  the  positive  X  axis  to  the  line  segmentOQ. 

We  suppose  P  to  be  the  location  of  a  2 -component  seismometer 
consisting  of  a  vertical  component  oriented  downward  along  Z,  and  a  hori¬ 
zontal  component  oriented  in  the  positive  X -direction. 

A.  FINITE  NUMBER  OF  DISCRETE  UNCORRELATED  POINT  SOURCES 

Suppose  we  have  a  finite  number  N  of  single -mode  Rayleigh 

wave  point  sources,  located  at  azimuth  0, ,  0„ . 0  .  Let  X  (t)  be 

12  N  n 

the  inline  horizontal  motion  at  P  due  to  the  n**1  source;  X  (t)  is  positive  in 

n 

the  direction  -0  ,  the  direction  of  propagation.  Let  v  (t)  be  the  vertical 

^  th  ^ 

motion  at  P  due  to  the  n  source. 


❖ 


We  make  the  following  basic  assumptions: 


(1)  For  each  n,  X^(t)  is  a  stationary  random  process. 


(2)  The  autocorrelation  functions  ^  (t)  are  all  the 

n  n 

same  except  for  multiplicative  constants. 


(3) 

(4) 


All  of  the  crosscorrelation  functions 
are  zero. 


There  exists  a  function  k(t)  with  / 
such  that  v  =  k  ©  X  for  all  n,  i.  e.  , 


n  n 


^  X  X  ,  m  ^  n , 
m  n 

|k(t)|dt  <  <*> 


v„(t) 


k(s)  X  (t-s)ds 
n 


In  this  report,  the  term  "Rayleigh  wave"  is  used  to  refer  to  all  types 
of  elliptically  polarized  surface  waves  and  is  meant  to  encompass  all 
surface  waves  except  Love  (SH)  waves. 
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Let  us  define 


h  (t)  =  -  cos  9  X  (t) 
n  n  n 


(3.  1) 


N 

h(t)  =  £  V0'  (3.2) 


v(t) 


vn(t) 


Clearly  is  the  contribution  to  the  output  of  our  horizontal 
seismometer  due  to  t  e  nth  source,  and  h(t)  and  v(t)  are  the  outputs  of  the 
hoiizontal  and  vertical  seismometers,  respectively. 


We  proceed  to  calculate  the  coherence 

2 


C  = 


1»J' 

*  Ku 

vv  hh 


between  v  and  h. 


Our  assumption  '3)  is  equivalent  to  another  assumption: 


(3a)  ^  (f)  =  0,  rn  ^  n 

m  n 

Our  assumption  (2)  is  equivalent  to  another  assumption: 

(2a)  There  exists  an  autopower  spectrum  $  ,  and  non-negative 

N  N 

constants  p(9  ),  n  =  1,  2,  .  .  .  ,  N,  with  £  p(6  )  =  1,  such  that 

n  n=l  n 

\  =  p(0n)  $  (f)»  n  =  1,  2,  .  .  .  ,  N  (3.  3) 

n  n 

At  this  point,  because  of  a  limiting  process  to  be  carried  out 
later,  it  is  desirable  to  introduce  one  additional  assumption: 


(3.4) 
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We  have 


=  $ 


h  h  ‘cos  9  ).  ,  cos 
n  n  n  n  n  n 


a  \  =  cos  9  $  =  cos  9  p  (9 

o  A.  n  a  X  n  n 


n  n 


Also,  by  assumption  (4)  and  formula  (3.  1)  we  get 


cos  9  v  k  @  h 
n  n  n 


Hence  by  (2.  7), 


-  cos  9  $  =  K4.  .  =  K  cos  9  p  (9  )  i 

n  v  h  h  h  n  n' 

n  n  n  n 


where 


00 

K(f)  =  J  k(t) 


-i2nft  , 
e  dt 


is  the  Fourier  transform  of  k(t).  That  is, 


♦  ,  =  -  K  cos  9  p  (9  )  $ 

V  h  nrx  n 


n  n 


From  assumption  (4)  and  (2.7)  we  get 


$v  V  $k©  X  ,  k@  X 
n  n  n  n 


=  |k|2  t,  ,  =  Ik |2P(0 


x  x 

n  n 


n 


Now,  l  om  (3.  5)  and  assumption  (3a)  we  get 


hh 


N  /  N  2  \ 

=  L  ‘h  h  =(E,  cos  6„p(9„>) 

n=l  n  n  \n=l  / 


(  3.  5) 

)• 


(3.6) 

(3.7) 


(3.8) 


(3.9) 

(3.  10) 
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Similarly  from  8)  and  {3.  9)  we  get 

$  ,  =(  2^  cos  9  p(9  )|K$ 

vh  \n=l  n  n  / 

and 

/  N  \  2  2 

}vv  =(Ep<V)|k|  }  *  Ik  I  * 

Therefore  KT 


I'/  I'S,  cos  e  p(9  )  K  I 

vn  n-i  n  n 


$  $ 

’  N 

- 

hh  vv 

£ 

cos^  9  p  (9  )  $ 

Ik  I2  $ 

n=l 

n  n 

. 

(3.  11) 


(3.  12) 


or,  finally, 


/  N  \2 

\  53  cos  9  p(9  )/ 
\n=l  n  ny 


E  cos  9  p{8  ) 
n=l  n  n 


(3.  13) 


We  call  p(8  )  the  (discrete)  azimuthal  power  distribution  function. 

n  - - * - - 

Thus  (3.  13)  gives  the  coherence  C  between  vertical  and  horizontal  for  N  discrete 

single -mode  Rayleigh  sources  with  azimuthal  power  distribution  function  p(9  ). 

n 

B.  CONTINUOUS  SOURCE 

It  is  of  interest  to  derive  a  formula,  similar  to  (3.  13),  giving 
the  coherence  for  the  case  of  uncorrelated  single-mode  Raylei  SL  noise 
arriving  from  the  continuous  range  of  azimuth  -n  <  9  <  tt  ,  with  a  continuous 
azimuthal  power  density  function  P(8). 

The  meaning  of  the  underlined  phrase  requires  some  clarification. 
Throughout  this  study,  noise  arriving  from  a  continuous  range  of  azimuths  will 
be  regarded  as  a  limiting  case  of  the  problem  we  have  already  considered 
involving  a  finite  number  N  of  uncorre.’ated  sources,  where  N  -»  ®  and  the 
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azimuthal  separation  between  sources  -*  0.  We  will  take  the  limit  in  such 
a  way  that  the  total  vertical  power  remains  constant. 


Thus,  suppose  P( 9 )  is  a  non-negative,  Riemann-integrable 
function  defined  on  |  -tt  ,  tt|  such  that 

/  P(9)  d  9  =  1 

-TT 

and  suppose  $  (f)  is  an  autopower  spectrum  such  that 


(3.  14) 


/  *(£) 


df  =  1 


We  call  P(9  )  a  continuous  azimuthal  power  distribution 
function.  Intuitively,  we  will  think  of  P(9)  as  defining  the  distribution  of 
source  power  with  regard  to  azimuth,  for  a  single -mode  Rayleigh  source 
distributed  over  the  whole  range  of  possible  azimuths. 

Identifying  the  endpoints  of  the  inteval  |-n ,  nj,  let  -TT  = 

9m  <  9,<  9  <  ...<  0  <9  =  TT  and  define 

N  1  2  N-l  n 


L 


3n+  1 


P ( 9 )  d9  ,  n  =  1,  2,  .  .  .  ,  N-l 


P(Qn>  =\ 


Jn 


rei 

/  P( 9 )  d0,  n  =  N 

-'0 


N 


Then  ^  p(9  )  =  1,  and  if  we  have  uncorrelated  Rayleigh 

n=l 


n 


sources  at  the  azimuths  9^,  .  .  .  ,  9^  with  inline  horizontal  autopower  spectra 
p(9^)  $ . .  p(9^)  $,  then  the  coherence  will  be,  by  (3.13), 


i. 
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V--  J; 


l 


■  0  % 


K 


— ,  cos  9  p(9  ) 
n=l  n  n 


r  • 


,  9 


n 


N  , 

]C  cos  9  >  (9  ) 
n=l  n  n 


(3.  16) 


N  ,  9 

E  Jt 

n=  1  n 


n+  1 


cos  9  P(  6  )  d0 
n 


II 

SMz 

f  9n+  1 

'e 

n 

2  , 
cos  d 

n 

P  (9)  d9 

If  we  let  N  increase  in  such  a  way  that 

min 

1  £  n^N 

1 9  -9 

1  n+  1 

easily  shown  that 

[N  r 

6n+  1 

t 

L  Jt 

|_n=  1 

)  cos 

n 

9P(9)  d9 

c  9  Q  c 

r  *  *  •  *  n 

N  r 

E  i 

9„+i 

cos 

2  9  P  ( 9)  c.  0 

n=l  0 

n 

r 

*-n 


cos  9  P(9)  d9 


/ 

-TT 


cos  9  P(B)  d9 


In  this  sense  we  say  that  the  formula 

,2 


/; 

c  =  --  n 

,n 


cos  9  P(  0  )  d  0 


cos  9 P( 0)  d0 


-n 

gives  the  coherence  between  vertical  and  horizontal  for  a  two-component 
seismometer  in  a  single -mode  Rayleigh  field  witn  continuous  azimuthal 
power  distribution  function  P(9). 


(3.  17) 
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C.  CALCULATION  OF  COHERENCE  FOR  CERTAIN  DISTRIBUTIONS  P(c) 

1.  C  =  1 

Having  obtained  formulas  (3.  131  and  (3.  17),  we  shall  now  use 
these  formulas  to  calculate  coherence  for  several  examples  of  azimuthal 
power  distribution  functions.  The  coherence  C  between  v  and  h  ranges 
from  0  (no  prediction  capability)  to  1  (perfect  predictability).  A  single 
discrete  source  will  yield  C=  1,  provided  its  azimuth  0  is  not  equal  to 

TT 

±  y  .  More  generally,  (see  formula  (3.13)),C=1  whenever  only  two  sources 
are  present  and  their  azimuths  satisfy  9  =  -9  9  ^  ±  .  (Power  may  be 

X  L  X  Ci 

distributed  arbitrarily  between  these  two  sources;  e.  g.  ,  one  source  may 
have  zero  power.)  It  may  be  shown  that  this  is  the  only  situation  for 
which  C  =  1 . 

2.  C  =  0 

The  other  extreme,  C=0,  can  occur  in  many  ways.  The  most 
important  cases  for  which  C=  0  are  the  following: 

(1)  Symmetry  About  -y 

Formula  (3.  17)  shows  that  if  P(  ^  -  9)  =  P  ( ^  +  9) 

for  all  9  ,  then  C  =  0. 

(2)  Periodicity 

If  there  exists  an  integer  S  s  2  such  that  P(9  +  -^)  =  P(9) 
for  all  9,  then  C  =  0.  This  follows  from  (3.7)  and  the  fact  that 

E  cos  (£*  ♦  ().  0 
k=l  '  b  ' 

for  all  positive  integers  S  ^  2  and  real  numbers  6.  An  important  example 
which  satisfies  both  (1)  and  (2)  is  the  isotropic  distribution: 

p(e>  =  l 
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However,  it  is  possible  for  the  coherence  to  be  quite  large, 
even  for  certain  ra°'er  widely  distributed  sources,  as  we  shall  illustrate 
with  the  three  samples  which  follow.  (See  Figure  1 1 1  - 1 . ) 

3.  Discrete  Point  Sources  of  Equal  Power  Located  at  Azimuths  U -p  and  p  +  p 


Consider  the  discrete  azimuthal  power  distribution  function 
p(M*P)  =  1/2,  p(pi+p)  =  1/2.  For  this  case,  the  coherence  computed  by 
(3.  13)  is 


.2 


P 


-J  cos  (|l-p)  +  ~  cos  (M+  p)  J 

1  2  ,  ,1  2, 

—  cos  ((a-P)  +  —  cos  (pi+P) 


2  cos  [a  cos  p 
1  +  cos  2  pi  cos2p 


(3.  18) 


A  tabulation  of  the  prediction  error  1 -C  for  several  values  of  pi  and  p 

a ,  P 

is  given  in  Table  III  —  1 . 

4.  Uniform  Distribution  From  (i-P  to  |i  +  p 


Let 


P(B)  = 


2p 


,pi-p£0spi+p 


(  0,  otherwise 


This  represents  a  uniform  distribution  of  power  coming  from 
a  wedge  of  azimuths  with  width  2p  and  center  angle  pi .  By  (3.  17)  the 
coherence  is 


H  +  P 


1  p 

2p  - p 


cos 


e  d  e 


M,  p 


1  r^+  P  2 

-  '  cos 


2p  J 


e  d  9 


M  -  P 


(3.19) 


sin 


(M-f  p)  -  sin(pi  -  p)l 


2p  P  + 


sin(2pt  +  2p)  -  sin(2pi  -  2p) 


,  2  .2 

4  cos  pi  sin  p 

p(2p  +  cos  2pi  sin  2p) 


A  tabulation  of  1-C^  p  f°r  certain  values  of  pi  and  p  is  given 


in  Table  III -2. 
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Table  III-l 

VALUES  OF  1-C  FOR  TWO  DISCRETE  SOURCES  OF  EQU^L  POWER 
LOCATED  AT  AZIMUTHS,  H-  p  AND  U  +  p.  (u.p  IN  DEGREES) 


ssa 

0 

10 

20 

30 

60 

90 

5  ;  175 

0.  000 

0.000 

0.  301 

0.  003 

0.  022 

1.  000 

10  ;170 

0.  000 

0.  001 

0.004 

0.  010 

0.  085 

1.  COO 

15  ;165 

0.  000 

0.  002 

0.  009 

0.  023 

0.  177 

1.  000 

20  ;1 60 

0.  000 

0.  004 

0.017 

0.  042 

0.  284 

1  ooo 

25  ;1 55 

0.  000 

0.007 

0.  028 

0.  068 

0.  ^95 

1.  000 

30  ;  1 50 

0.  000 

0.  010 

0.  042 

0.  100 

0  500 

1. 000 

35  ;1 45 

0.  ooo 

0.015 

0.061 

0.  140 

0.  595 

1. 000 

40  ;  1 40 

0.  009 

0.  021 

0.085 

0.  190 

0.  678 

1. 000 

45  ;1 35 

0.  000 

0.  030 

0.  117 

0.250 

0.750 

1. 000 

50;  130 

0.  000 

0.  042 

0.  158 

0.  321 

0.  810 

1.000 

55 ; 125 

0.  000 

0.  060 

0.  213 

0.  405 

0.  860 

1.  000 

60 ; 120 

0.  000 

0.  085 

0.284 

0.  500 

0.  900 

1. 000 

65  ;  1 1  5 

0.  000 

0.  125 

0.  379 

0.  605 

0.  932 

1.  000 

70  ;  1 1 0 

0.  000 

0.  190 

0.  500 

0,  716 

0.  958 

1. 000 

7 5  ;  105 

0.  000 

0.  302 

0.  649 

0.  823 

0.  977 

1.  cco 

80  ;  100 

0.  000 

0.  500 

0.  810 

0.  914 

0.  990 

1.  009 

85  ;95 

0.000 

0.  802 

0.  945 

0.  977 

0.  997 

1. 000 

90 

1, 000 

1.  000 

1. 000 

1.  000 

i.  ooo 
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VALUES  OF  1 -C  FOR  UNIFORM  DISTRIBUTION 
FROM  U-p  TO  M+p  (n.p  IN  DEGREES) 


kS 

0 

10 

20 

30 

60 

10 

5 

0.  000 

0.  000 

0.  000 

0.  001 

0.  008 

1.  GOO 

15 

0.  000 

0.  001 

0.  003 

0.  008 

0.  065 

1. 000 

25 

C.  001 

0.003 

0.  009 

0.  022 

0.  164 

1.  000 

35 

0.  003 

0.  007 

0.  020 

0.045 

0.  284 

1. 000 

45 

0.  009 

0.  016 

0.  038 

0.  078 

0.  405 

1. 000 

55 

0.  022 

0.  033 

0.065 

0.  122 

0.  518 

1. 000 

65 

0.  046 

0.  060 

0.  104 

0.  181 

0.  616 

1. 000 

75 

0.  086 

0.  105 

0.  161 

0.  254 

0.  699 

1.  000 

85 

0.  148 

0.  171 

0.  238 

0.  343 

0.  768 

1.  000 

95 

0.  238 

0.  264 

0.  336 

0.  444 

0.  824 

1. 000 

105 

0.  357 

0.  382 

0.452 

0.  553 

0.  870 

1.  000 

115 

0.  496 

0.  518 

0.  578 

0.  662 

0.  907 

1. 000 

125 

0.  641 

0.  657 

0.  702 

0.  763 

0.  936 

1.  000 

135 

0.  771 

0.  782 

0.810 

0.  849 

0.  959 

1.  000 

145 

0.  874 

0.  879 

0.  894 

0.  915 

0.  976 

1.  000 

155 

0.  943 

0.  945 

0.  952 

0.  961 

0.  989 

1.  000 

165 

0.  982 

0.  983 

0.  985 

0.  987 

0.  996 

1.  000 

175 

0.  998 

0.  998 

0.  999 

0.  999 

1.  000 

1.  000 

180 

1. 000 

1.000 

1.  000 

1.  000 

1. 000 

1.  000 
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5.  Cos  Distribution 


Consider  the  continuous  azimuthal  power  density  function 


P(G)  = 


(  1  2  /tt 

)  p  \2p 

v  0,  otherwise 


(6  -  MM,  U  -  P  s  6  *  U  +  P 


Using  formula  (3.  17),  we  find.,  for  p  ^  —  ,  p  ^  n 
/  2  n2  \  2  .  2 

— Tj  cos  ^  sin  P 


2P  P 


[  n_ 
P  +  — - 


L  2n2-3p2 


cos  2 (i  sin  2p 


(3.20) 


n 


2 

cos 

1+2  cos2  |i 


(3.21) 


2 

cos  pi 


A  tabulation  of  1-C^  ^  for  several  values  of  M,  p  is  given 


in  Table  III -3. 
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Table  III -3 


VALUES  OF  1-C  FOR  COSZ  DISTRIBUTION 
(b,  P  IN  DEGREES) 


B 

0 

10 

20 

30 

60 

90 

5 

0.  000 

0.  000 

0.  000 

0.  000 

0.  003 

1.  000 

15 

0.  000 

0.  000 

0. 001 

0.  003 

0.  026 

1.  000 

25 

0.  000 

0.  001 

0.  004 

0.  008 

0.  070 

1.  000 

35 

0.  001 

0.  002 

0.  00 1 

0.  017 

0.  129 

1. 000 

45 

0.  002 

0.  005 

0.  013 

0.  029 

0.  199 

1.  OwO 

55 

0.  005 

0.  009 

0.  021 

0.  044 

0.  273 

1. 000 

65 

0.  010 

0.  016 

0.  0  32 

0.  °64 

0.  348 

1.  000 

75 

0.  019 

0.  026 

0.  048 

0.  089 

0.  420 

1.  000 

85 

0.  031 

0.  040 

0.  068 

0.  1 18 

0.  488 

1. 000 

90 

0.  039 

0.  049 

0.  080 

0.  135 

0.  520 

95 

0.  049 

0.  060 

0.  094 

0.  154 

0.  550 

1.  000 

105 

0.  073 

0.  086 

0.  126 

0.  195 

0.  607 

1.  000 

115 

0.  106 

0.  120 

0.  165 

0.  24i 

0.  657 

1. 000 

125 

0.  146 

0.  162 

0.  211 

0.  293 

0.  703 

1 . 000 

135 

0.  196 

0.  213 

0.  265 

0.  349 

0.  743 

1.  000 

145 

0.  253 

0.  271 

0.  324 

0.  409 

0.  779 

1. 000 

155 

0.  319 

0.  336 

0.  389 

0.  471 

0.  811 

1.  000 

165 

0.  389 

0.  406 

0.  456 

0.  534 

0.  839 

1. 000 

175 

0.  463 

0.  479 

0.  525 

0.  595 

0.  864 

1. 000 

180 

0.  500 

0.  515 

0.  558 

0.  625 

0.  875 

1.  0.^0 
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SECTION  IV 


I 


I 

i 
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EFFECT  OF  ADDITIONAL  RANDOM  UNCORRELATED  NOISE 


The  coherence  between  vertical  and  horizontal  for  the  case 
we  have  considered  in  Section  III  will  be  reduced  if  the  seismometer  outputs 
v  and/or  h  contain  an  additional  component  of  random,  uncorrelated  noise; 
e.g  .  ,  amplifier  noise. 

Consider  the  following  quite  general  problem:  we  are  given 
two  stationary  random  processes  v  and  h,  with  coherence  C,  v^  =  v  +  z, 
h  =  h  +  x,  where  z  and  x  are  stationary  random  processes  uncorrelated 

to  v,  h,  and  uncorrelated  to  each  other,  then  what  is  the  coherence  C* 

1  ,  ,  1  „ 
between  v  and  h  ? 


We  have 


$ 


1  1 
v  h 


$ 


vh 


v  v 


vv  zz 


(4.1) 


We  see  from  equation  (4.  2)  that  the  effect  on  coherence  depends  upon  the 

ratios  $  /$  and  $  /$.  .  .  In  the  case  that  v  and  h  are  the  vertical 

zz !  vv  xx/  hh 

or  horizontal  components  of  displacement  due  to  a  single-mode  Rayleigh 

field  with  some  azimuthal  power  distribution,  P(9),  and  if  we  conceive 

that  $  ,  $  ,  and  $  remain  fixed  as  we  vary  P(8),  we  see  that  the 

zz  vv  xx 
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greatest  reduction  in  coherence  results  when  P(9)  is  strongly  concentrated 


TT 


near  ±  —  ,  for  then,  by  equation  (3.  10),  $  will  be  very  small.  For  example, 
2  hh 

if  <  $  ,  then  C1  "  1/2  C,  by  (4.  2). 

hh  xx 


true  that 


For  any  real  numbers  CL,  3  with  0  -  0L  £  1,  0  s  9  £  1 ,  it  is 

(r H)  (ih) 


>  l  -  a  -  3 


Therefore,  if 


$  s  $  $  s  $ 

zz  vv  xx  hh 


then 


1  - 


zz 


vv 


$ 

XX 

r 

$ 

hh  J 

$  and  $ 

zz 

TT 

(4.3) 


(4.4) 


XX 


(4.  3)  shows  that  if  on  each  channel  the  power  of  amplifier  loise  is  less  than 
1  percent  of  the  power  of  seismic  noise,  then  the  coherence  is  reduced  by 
less  than  2  percent. 


We  may  conclude  that  for  the  case  of  single-mode  Rayleigh 
waves  with  an  azimuthal  power  density  P(9),  the  coherence  between  v  and 
h  will  not  be  seriously  affected  by  the  addition  of  a  realistic  amount  of 
amplifier  noise  on  both  channels,  unless  P(9)  is  concentrated  near  ±  — -  . 
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SECTION  V 


LOVE  WAVES 


We  shall  discuss  in  this  section  the  limitation  upon  prediction 
capability  imposed  by  the  addition  of  Love  (SH)  waves  on  the  horizontal 
channel.  It  will  always  be  assumed  that  the  Love  waves  are  uncorrelated 
with  the  Rayleigh  waves. 


A.  GENERAL  FORMULA  FOR  EFFECT  OF  LOVE  WAVES  ON  COHERENCE 

Suppose  we  have  N  sources,  at  azimuths  9^ . 9  ,  each 

source  generating  Rayleigh  noise,  with  components  v  ,  >.  ,  h  as  in  Section 

th  n  n  n 

III.  Suppose  also  that  the  n  source  generates  Love  wave  noise  with  motion 

i  in  the  direction  9  +  tt  /  2  . 

n  n 

In  addition  to  the  assumptions  of  Section  III,  paragraph  A, 

we  assume 


and 


(1)  For  each  n,  jn(t)  is  a  stationary  random  process 

(2)  There  exists  an  autopower  spectrum  Q(f)  and  constants 

N 

q(9  )  such  that  ^  q(9  )  =  1  and  $.  .  =  q(9  )  Q, 

n  n=  1  n  J  J  n 

n  1  n  n 

n  =  .! ,  2 . N. 


(3)  The  crosspower  spectra  $.  =  0,  m  ^  n  and 

J  J 
m  n 

The  crosspower  spectra  $.  ^  =  0  for  all  m,  n. 

^m  n 


Let  us  define 


w  (t)  =  sin  9  ;  (t) 

n  n  Jn 


(5.  1) 


and 


N 

w(t)  =  £  w  (0  (5.2) 

n=  1  n 


V-l 


Clearly  w  is  the  contribution  to  the  output  of  the  horizontal 
^  th 

seismometer  due  to  Love  waves  from  the  n  source,  and  w(t)  is  the  total 
Love  component  on  the  horizontal. 

Let  L(f)  be  defined  by  the  relationship 


n  (f)  =  L(f)  $  (f) 


(5.  3) 


Then  (cf. 

$ 

vvw 


3.  10)  we  find 

£  .  2  . 
=  >  sm  a 

n=  1 


n 


q(R  )  n 

n 


(5.4) 


2 

=  ^  sin  0  q  ( 9  )  •  L  •  $ 
n=  1  n  n 


by  Equations  (3.  10)  -  (3.  12),  (5.  4).  Hence 

N  l2 

/-j.  cos  9  p(9  ) 

„  n=  1  nr  n  _ 

L  =  N  2  M 

iLd  cos  9  p(9  )  +  L  £  sin2  9  q(0  ) 
n= 1  n  n  n=l  n  n 

We  call  q(9  )  the  (discrete)  azimuthal  power  distribution  function 
n 

for  the  Love  waves.  As  before,  we  can  obtain  a  formula  similar  to  (5.  6) 
giving  the  coherence  for  a  continuous  azimuthal  power  distribution  Q( 9 ) 
for  the  Love  waves: 


TT 

L 


cos  9  P( 6  )  d9 


L 


cos  ~P(~)d9  +  L'l  _  sin^0Q(9)d9 

TT 


(5.7) 


Note  that  (5.  6)  and  (5.  7)  reduce  to  (3.  1  3)  and  (3.  17)  respectively, 
on  setting  L  =  0. 

B.  ISOTROPIC  LOVE  WAVES 

It  is  best  to  regard  the  case  of  isotropic  Love  waves 
(ate) »  i)as  a  special  case  of  the  situation  discussed  in  Section  IV.  That  is, 
the  Love  wave  energy  on  the  horizontal  is  simply  uncorrelated  additive  noise 
on  the  predicting  channel,  and  its  effect  upon  coherence  is  mathematically 
indistinguishable  from  that  of,  say,  amplifier  noise. 


(4.  2)  that 


Thus,  if  C  is  the  coherence  between  v  and  h,  we  have  by 


hh 


$  +  $ 

hh  ww 


(5.8) 


If  §  is  large,  there  is  a  substantial  reduction  in  coherence, 
ww 


For,  recall  (cf.  3.  10)  that 


hh 


TT 


J  cos  0  P(6)  d9 

-TT 


$  £  $ 


Also,  by  (5.4) 


Therefore 


ww 


/. 


sin  9  •  —  d9 
2rr 


=  1  Li 
2 


hh 


$ 


hh  +  $ 


ww 


i  +  ~  Li 

Cm 


2  +  L 
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Hence, 


isotropic  Love  waves 


(5.9) 


Formula  (5.  9)  implies,  for  example,  that  if  L  =  2,  then  the 

coherence  C  is  never  greater  than  1/2. 

.Li 

C.  LOVE  WAVES  WITH  SAME  DISTRIBUTION  AS  THE  RAYLEiGH  WAVES 

If  the  Love  wave  field  is  non -isotropic,  formula  5.  8  still 

holds.  But  in  this  case  it  may  be  misleading  to  regard  the  Love  waves 

simply  as  additive  uncorrelated  noise,  because  $  depends  upon  the 

ww 

orientation  of  our  horizontal  seismometer.  That  is,  by  pointing  our 

horizontal  seismometer  in  the  proper  direction,  it  may  be  possible  to 

make  $  quite  small, 
ww 

To  illustrate  this  phenomenon,  consider  the  situation  in 
which  the  Love  waves  have  the  same  azimuthal  power  distribution  as  the 
Rayleigh  waves;  i.  e.  , 


Q(9)  =  P(  9 ) 
Then  (5.7)  becomes 


TT 


f 

■■'-TT  cos  9  P(9)  d  9 


(5.11) 


TT 


-TT 


2  rn 

tos  9P(8)d0  +  L  I  sin*  9p(0)d8 

-TT 


If  in  the  distributions  P(9)  considered  in  Section  III,  paragraph  C, 
we  specify  the  values  of  p  and  L,  then  the  coherenct  CT  calculated  by  (5.  11)  is 
maximized  for  p  =  0.  For  pi  =  0,  we  get  surprisingly  higi.  coherence,  even 
for  L  as  large  as  10  (See  Tables  V-l,  V-2).  For  example,  for  a  iform 
dist  .bution  from  -25°  to  +25°  and  L  =  10,  the  prediction  error  is  oniy 
0.  395  (Table  V-2). 


-4 


Table  V-l 


VALUES  OF  1  -C  FOR  DISCRETE  SOURCES  AT  AZIMUTHS  ±  P 
WITH  INLINE  LOVE  POWER  =  L  x  INLINE  RAYLEIGH  POWER 


(Di-Arif's) 

L  =  0 

L  =  0.  1 

L  =  1 . 0 

— 

L  *  10.  0 

1  7  5 

0,  000 

0.  001 

0.  008 

0.  07  1 

1  =>,  165 

0.  000 

0.  007 

0.  067 

0.  418 

25,  155 

0.  000 

0.  0i! 

0.  179 

0.  685 

15  145 

0.  000 

0.047 

0.  529 

0.  8  51 

45;  1  55 

0.  000 

0.  091 

0.  500 

0.  909 

55;  1  2  5 

0.000 

0.  169 

0.  671 

0.  955 

65;  115 

0.  000 

0.  3 1  5 

0.  821 

0.  979 

7  5;  105 

0.  000 

0.  582 

0.  933 

0.  993 

85;  95 

0.  000 

0.  929 

0.  992 

0.  999 

Table  V-2 

VALUES  OF  1-C  FOR  UNIFORM  RAYLEIGH  DISTRIBUTION  FROM  -p  TO  +  p 
WITH  INLINE  LOVE  POWER  =  L  x  INLINE  RAYLEIGH  POWER 


0 (Degrees) 

o 

ii 

L  =  0.  1 

o 

If 

L  =  10.  0 

5 

0.  000 

0.000 

0.  00  : 

0.  025 

15 

o.  ooo 

0.  002 

0.  023 

0.  187 

25 

0.  001 

0.  007 

. .  ut><: 

0.  395 

35 

0.  00  3 

0.016 

0.  118 

0.  568 

45 

0.  009 

0.  0  30 

0.  189 

0.  692 

55 

0.  022 

0.  054 

0.  271 

0.  779 

65 

0.  046 

0.091 

0.  361 

0.  840 

75 

0.  086 

0.  143 

0.455 

0.  882 

85 

0.  148 

0.217 

0.  549 

0.  914 

95 

0.  238 

0.  ‘.14 

0.  639 

0.  937 

105 

0.  357 

0.431 

0.  722 

0.  955 

115 

0.  496 

0.  5c  1 

0.  796 

0.  968 

125 

0.  641 

0.  689 

0.  859 

0.  978 

1  35 

0.  77! 

0.  802 

0.  910 

0.  986 

145 

0.  874 

0.  890 

0.  949 

0.  99£ 

155 

0.  943 

0.  950 

0.  976 

0.  996 

165 

0.  982 

0.  984 

0.  992 

0.  999 

175 

0.  998 

0.  998 

0.999 

1. 000 

180 

t  . 

1. 000 

1.000 

1. 000 

1 . 000 
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SECTION  VI 


EFFECT  OF  GAIN  FLUCTUATION 

In  the  preceding  sections  we  have  assumed  implicitly  that 
the  amplitude  gains  of  our  two  seismometer  s  are  constant  functions  of  time. 
In  practice,  however,  these  gains  vary  slowly  with  time,  and  are  never 
known  precisely.  This  results  in  a  lower  observable  coherence  between 
the  two  channels  than  is  indicated  by  our  previous  equations.  Fortunately, 
this  effect  is  not  serious  if  the  gain  fluctuations  are  reasonably  small  and 
veil  -behaved. 

For,  let  r  =  r(t)  a.'/!  s  -  sft)  be  the  gains  of  the  veit‘>'l 
and  horizontal  seismometers  respectively.  We  assume  that  ;  and  s  are 
non-negative  variables  satisfying: 

(1)  r(t)  and  s(t)  are  statistically  independent  stationary 

random  processes,  with  first  order  probability 

functions  Q  and  Q  ,  respectively 
r  s 

(2)  During  any  short  interval  of  time  during  which 
we  continuously  record  data  for  use  in  an  experi¬ 
mental  determination  of  cp  ,  cp  ,  ,  or  cp.  ,  ,  the 

vv  vh  hh 

variations  in  the  gains  are  so  slight  that  we  may 
consider  the  gains  to  be  constant  throughout 
t'-zt  interval 

With  these  assumptions,  suppose  now  that  a  certain  noise 

situation  persists  for  many  days;  that  is,  the  actual  correlations  cp  , 

vv 

cp^,  and  vp  ,  and  the  azimuthal  power  distribution  function  P(9)  remain 

unchanged  over  an  extended  period  of  time.  During  this  period  let  us 

repeatedly  use  our  two  seismometers  to  measure  cp  ,  cp  ,  and  cp,  ,  .  The 
'  vv  vh  hh 

gains  r  and  s  are  varying  dur;ng  this  period,  but  in  accordance  with 
assumption  (2),  we  make  the  approximation  that  r  and  s  remain  constant 


VI-1 


throughout  each  interval  of  measurement.  Then,  the  observed  correlations 
determined  by  measurement  beginning  at  time  t  are 


2  2 
r(0  ^vv’  ^vh’  and  ^hh 


Taking  the  averages  of  these  correlations  for  many  different 

values  of  t,  we  obtain  the  average  observed  correlations  <b  ,  <t>  .  ,  and 

- ° -  vv  vh 

*hh-  Thus 


$  =  r(t)  fip 

vv  vv 


(6.1) 


*vh  -  r!t!  -W  '•’vh 


+hh  =  s(,)  'Ohh 


H 

We  define  C  , 


the  observed  coherence,  by 

,2 

i 


(  r (t)  a (t) )  |  *vh| 


r,t>2  8(t'2  »vAh 


(6.  2) 


Therefore, 

c’ 

Ht)7 

where  C  is  the  coherence  between  v  and  h  that  would  be  observed  by 
seismometers  having  no  gain  fluctuation. 


(6.  3) 


The  mathematical  significance  of  C11  lies  in  the  fact  that 
1  -C"  is  the  relative  prediction  error  corresponding  to  the  optimum  Wiener 
prediction  filter  transfer  function  G(f)  to  be  applied  to  h  for  predicting  v  in 
the  variable  gain  situation. 
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In  greater  detail:  let  g(t)  be  a  convolution  operator,  G(f)  = 
Fourier  transform  of  g,  and  define  e  =  v  -  (g©h).  Let  <J>  be  the  average 


ee 


observed  autocorrelation  function  of  e  (cf.  the  definition  of  <b  ,6  and  6  ) 

vv  Wh  rhh 

and  let  P  (f)  be  the  Fourier  transform  of  <b 

e  Tee 

A(f) 

Take  G(f)  =  -  - : "  ,  where  A  and  B  are  the  Fourier  transforms  of 

£5(1) 

4>vh  and  4>hh  respectively;  it  may  be  shown  that  this  choice  of  the  filter 
transfer  function  G  is  optimal  in  the  sense  that  it  minimizes  the  average 
error  power  m 

f  P  (f)df 

J  ~  e 


Furthermore,  we  find  that  for  this  choice  of  G, 


00 

/  P  (f)  df  =  D(f)  [l-  c"(f)  ] 

—  00  c  I  J 


where  D(f)  is  the  Fourier  transform  of  <j> 


vv 


Now,  by  ergodicity  of  the  time  series  r(t)  and  s(t),  (6.  3). 


may  be  written 


.// 


Iff 


Q  (x)  Q  (y)  dx  dy 
r  s 


r®  2  r00  2 

/  Q  (x)  dx  /  Q  (y)  dy 

Jo  r  -b  s 


(6.  4) 


or,  since  r(t)  and  s(t)  are  statistically  independent, 

|2  r  . ®  ,2 


„  ll  Qr,x)dx]  [{  Q°(y,dy 


C.  = 


/  Qr  (x)2  dx  •  J  Qg(y)2  dy 


•  C 


(6.5) 


2  2/12  2 

^  a  /  V  +  ° 

r  r '  '  s  s 


where  |i  and  |i  are  the  means  of  the  distributions  Q  and  Q  respectively, 
r  s  r  s  7 

and  and  are  the  corresponding  variances. 
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it  is  clear  from  formula  (6.  5)  that  the  effect  on  coherence 

2  2 

due  to  gain  fluctuation  is  small  if  the  variances  o  and  a  are  small  com- 
2  2  r  s 

pared  to  and  To  illustrate  the  quantitative  effect  of  gain  fluctuation 

on  coherence,  we  consider  two  examples  of  gain  probability  distributions. 
Example  1.  Uniform  Distribution  from  0  to  R 

Let  R  b*»  an  arbitrarily  large  positive  number  and  let 


Q  (x)  =  Q  (x) 
r  s 


]_ 

R,  0  s  x  s  R 
0,  otherwise 


Then  compute  from  (6.  5)  that 


C 


a 


C 


Thus  in  the  rather  extreme  case  that  the  gains  are  evenly  distributed 
between  0  and  some  large  number,  the  coherence  is  reduced  by  less  than 
a  factor  of  2. 


Example  2.  Log  (gain)  Normally  Distributed 


Suppose  that  the  random  variables  log  r  and  log  s  both  are 

2 

norma:  i  distributed  with  mean  0  and  variance  O  .  Then 


M 

l 


Also, 


2  2 

|i  +  o 
r  r 
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Therefore,  by  (6,  5) 


(6.6) 


values  of  0 


2 


In  Table  VI  -1, 


and  e 


are  tabulated  for  several 


We  see  from  Table  VI -1  that  in  Example  2  (which  is  probably 

a  more  realistic  gain  distribution  than  Example  1)  a  severe  loss  in  coherence 
2 

results  if  0^  is  allowed  to  be  as  large  as,  say,  1.  It  should  be  possible, 
however,  through  a  program  of  peridical  calibration,  to  insure  that  the 
seismometer  gain  fluctuations  remain  within  acceptable  limits.  This  will 
be  necessary  if  one  is  to  implement  a  processor  in  which  the  prediction 
filter  applied  to  h  for  predicting  v  remains  unchanged  for  an  extended 
period  of  time. 


Table  VI -1 

EFFECT  OF  GAIN  FLUCTUATION  ON  COHERENCE  FOR  log  fgain) 
NORMALLY  DISTRIBUTED,  MEAN  =  0,  VARIANCE  =  a2 


a2 

r 

2 

a 

r 

-20  2 
e 

0 

1 

0 

1 

0.01 

1.0051 

0.0101 

0.9802 

0.1 

1.0513 

0. 1162 

0.8187 

0.3 

1.1618 

0.4722 

0.  5488 

0.5 

1.2840 

i.  0696 

0.  3679 

1.0 

1.6487 

4. 6708 

0. 1353 
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SECTION  VII 


MULTIMODE  NOISE 


In  Sections  III  through  VI  we  have  discussed  the  coherence 
between  vertical  and  horizont?!  seismometers  at  point  location  in  a  single¬ 
mode  noise  field,  with  the  possible  addition  of  Lrve  waves  on  the  horizontal. 
For  the  same  two-channel  system,  we  shall  now  obtain  formulas  for  coherence 
between  v  and  h  in  the  more  general  case  that  not  one,  but  several  modes  of 
noise  are  present. 


As  in  Section  III,  assume  we  have  a  two-component  seis¬ 
mometer  located  at  the  origin  O  of  an  X-Y  coordinate  system  in  the  plane. 

Let  9^.  . 6^  be  az*mut*ls;  for  each  n  =  1,  2,  .  .  . ,  N,  we  have  a 

noise  generator  at  azimuth  9  . 


Let  v  and  X  be  the  respective  outputs, due  to  noise  from 
th  n  .  n 

the  n  source,  of  a  vertical  seismometer  at  P  and  an  inline  horizontal 
seismometer  at  P. 


We  assume  that 


v 

n 


M 


■  E 


m 

v 

n 


X 


m 


n 


(7.1) 


where  v  and  X  are  respectively  the  vertical  and  inline  horizontal 
n  n 

components  of  m**1  noise  mode  from  the  n  source,  m  =  0,  1,  2,  .  .  . 


Also,  we  assume  that  the  n  source  generates  Love  (SH) 

noise  with  motion  j  in  the  direction  0  +  . 

n  n  2 


M. 


We  make  the  following  basic  assumptions  about  vm ,  Xm 

n  n 

and  j  : 


Notation:  numerical  subscripts  refer  to  the  source;  superscripts  refer 
to  the  mode. I 
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(1)  i~nd  j  are  stationary  random  processes,  for  all  n,  m. 

(2)  For  each  m  =  0,  1,  2,  .  ,  ,,  M,  there  exists  an  autopower 

spectrum  $m,  and  constants  pm(9  ),  such  thct 

n 

N  _ 

Epm(en)  =  i 

n=l 


’  =  Pm(6  )  ,  n  =  1,  2 . N 

Xm  Xm  n 

n  n 


There  exists  an  autopower  spectrum  0  and  constants 


q(0  )  such  that 
n 


N 

*  1 


$i  i  =  q(8n}  °*  n  =  1.  2.  •  •  ..  N 
JnJn 

Xnm  and  X^V  are  uncorrelated  unless  n  s  v  and  m  =  H. 


X  and  j  are  uncorrelated  for  all  n,  m,  p.  . 
m  p 


For  each  m,  there  exists  a  function  k  (t)  with 


r  ikm  (t)  i  dt 


<  CD 


such  that 


m  _  m  /jN  « m  .  _  _  _  _ 

v  =  k  ©  X  for  all  n  =  1,  2,  ,  ,  ,  ,  N 

n  n 
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Having  made  these  assumptions,  let  us  define 


h  (t)  =  -cos  6  >.  (t) 

n  n  n 


w  (t)  =  sin  9  j  (t) 


N 

h(t)  =  V  h  (t) 
n 

N 

w(t)  =  V*  w  (t) 
",  n 
n=l 

N 

v(t)  =2  v  (t> 


n=  1 


n 


(7.2) 


Clearly  is  the  contribution  to  the  output  of  our 

th 

horizontal  seismometer  due  to  the  n  source  (w  is  the  Love  wave  contribution); 
h(t)  +  w(t)  and  v(t)  are  respectively,  the  total  horizontal  and  vertical  seismo¬ 
grams. 

We  proceed  to  find  the  c..  orence 


C 


I  h  + w  I 

$  K 

vv  h  +  w,  hlw 


between  v  and  h  +  w. 

As  before,  let  K  be  the  Fourier  transform  of  km,  m  =  0, 
1,2,  .  .  .,  M,  and  let  L(f)  be  defined  by 


M 

0(f)  =  L(f)S  *m(f) 

m=0 

|l-.(f)  is  thus  a  measure  of  the  ratio  of  inline  Love  power  to  total  inline 
SV+P  power  . 


(7.4) 
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Omitting  the  calculations,  which  are  but  a  slight  generalization  of 
the  corresponding  calculations  in  Sections  III  and  V,  we  find 


and 


■£.  t 


*  .  ,  -E  D  cos  9  pm(0  )  Km  $m 

v,h+w  ^1  m=0  n  n 


vv 


-t  I 

m=0 


Km|2im 


N  M 

^h+w,  h+w  ", 

n=l  m=0 


cos^  9  pm  (9  )  $m 

n  n 


rs  2 

+  V*  sin  9  q(9  )  Q 
n=  1  n  n 


(7.5) 


One  can  now  substitute  (7.  5)  into  (7.  3)  and  obtain  the  general 
multimode  coherence.  However,  we  shall  now  restrict  ourselves  to  a  simpler 
case;  namely,  we  make  one  additional  assumption: 


(5)  All  of  the  azimuthal  power  distribution  functions  are 

,  .  0  1  m 

equal,  i,  e.  ,  p  =  p  =  .  .  ,  =  p  =  q. 


Letting  p  be  the  azimuthal  power  distribution  function  common 
to  all  modes,  in  accordance  with  assumption  (5),  we  find  from  (7.  3),  (7,4)  and 
(7,5)  that  the  coherence  is: 


C  = 


/  n  y 

(  E. cos  9  p(0  )  ) 
\n=l  n  n  / 


(7.6) 


^  T-m  *m 


m=0 


K  r 


N 

E 


N 


2  M 

— .  cos  9  P(9J  +  L  E  sin  9  p(9  )  E  $m*Eft 
n=  1  n  n  n=  1  n  n  m=0  m=0 


K 


m 


,m 


VII -4 


As  before,  we  can  write  a  similar  formula  for  the  case  of  a 

2  (7.7) 


continuous  azimuthal  power  distribution  function  P(^); 


C  = 


r 

/  cos  *5  P{ 9)  d  9 


M 

Jo 


£„  Km  im 


_TT 


TT 

f  2 

M 

M 

2 

[  sin  9  P(9)  d9 

£ 

§m  •  E 

Km 

-TT 

o 

ii 

£ 

o 

II 

£ 

Formula  for  multimode  coherence,  all  modes  having  same 
azimuthal  power  distribution  function  P(6). 

Comparison  of  (7.7)  with  (5.  11)  shows  that  (Multimode 
coherence)  =  Q  (f)  X  (Single -mode  coherence),  where 


Q 


M 

n5o 


Km  $m 


M 

£  * 

m=0 


m 


M 

•I0  ik 

m=0 


m 


tm 


(7.8) 


is  a  function  depending  only  on  the  autopower  spectra  and  modal  horizontal- 
vertical  transfer  functions  K™.  Since  we  have  already  investigated  single -mode 
coherence  in  Sections  III  through  VI,  the  study  of  multimode  coherence  for  the 
case  that  all  modes  he  e  the  same  azimuthal  power  distribution  function  reduces 
to  a  consideration  of  the  function  Q(f). 


Detailed  examination  of  the  behavior  of  Q(f)  for  realistic 
choices  cf  the  $m's  and  Km' s  has  not  been  carried  out.  However,  we  can 
make  a  few  elementary  remarks. 

First,  it  should  be  noted  that  the  presence  of  more  than  one 
mode  always  results  in  a  reduction  in  coherence  from  the  single-mode  case. 

For,  it  follows  from  Holder's  inequality  that  (for  real  frequencies)  0  £  Q  (f)  £  1  „ 
Furthermore,  assuming  all  $  to  be  non-zero,  we  have  Q(f)  <  1  unless  K  =  K1  * 
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In  a  physical  situation,  it  is  extremely  unlikely  that  all  the 

transfer  functions  Km  would  be  equal,  or  nearly  so.  Approximately,  we 

might  expect  that  there  would  exist  positive  constants  R  ,  R  ,  .  .  R 

0  1  M 

such  that  over  a  fairly  broad  band  of  frequencies  we  have 

Km  =  ±iR  ,m  =  0,  1,...,M 
m 

Take  as  an  example  the  case  M  =  1  (two  modes),  and  let  =  $  .  Let  R 
be  any  positive  number.  Table  VII-1  shows  values  of  Q  computed  by 
(7.8)  for  several  choices  of  K^,  K^. 


Table  VII-1 

VALUES  OF  Q  FOR  TWO-MODE  CASE,  $°  =  $ 1 


Ko 

Ki 

Q 

+  i  R 

+  i  R 

1 

+  i  R 

+  2  i  R 

9/10 

+  i  R 

+  5  i  R 

9/13 

+  i  R 

+  b  i  R 

(b+l)2/2(b2  +  l),  b>0 

+  i  R 

-  i  R 

0 

Table  VII-1  illustrates  a  general  principle,  which  is 
roughly  as  follows:  the  reduction  in  coherence  due  to  additional  modes  is  most 
severe  in  the  case  that  some  of  the  Km’s  represent  +90°  phase  shift  and  others, 
a  -90°  phase  shift;  for  frequencies  where  all  Km,s  have  the  same  phase  response, 
the  effect  is  less  marked.  Indeed,  for  the  two-mode  case,  we  get  as  low  as 
Q  s  0  if  and  K*  have  opposite  phase  responses,  whereas  Q  if  and  K* 
have  the  same  phase  responses. 
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Given  an  azimuthal  power  distribution  function  P(9)  it  is 
apparent  that  if  several  modes  are  present  with  approximately  equal  power 
in  each  mode,  then  the  coherence  is  likely  to  be  substantially  reduced  from 
that  for  a  single  mode  with  the  same  distribution  P(^).  On  the  other  hand, 
in  closing  this  section  it  should  be  remarked  that  it  also  follows  from  (7.8)  that 
if  one  mode--s'.y,  the  fundamental  Rayleigh  mode  - -accounts  for  all  but  a  small 
amount  of  the  total  power,  then  there  will  be  very  little  reduction  in  coherence 
from  the  single-mode  case. 


VII -7/ 8 


BUNK  PAGE 


SECTION  VIII 
DISCUSSION 


In  the  preceding  sections,  mathematical  formulae  have  been 
derived  for  the  noise -rejection  properties  of  a  2-component  seismometer, 
consisting  of  1  vertical -component  and  1  horizontal-component  instrument 
at  the  same  point.  For  most  applications,  this  system  may  be  considered 
equivalent  to  a  3-component  seismometer.  It  is  possible  to  design  a 
processing  system  wherein  different  filters  are  applied  to  the  outputs  of 
2  horizontals  and  the  noise  estimate  is  obtained  by  summing  the  filter 
outputs.  This  case  is  not  covered  by  the  development  of  Part  I,  but  the 
simpler  case,  in  which  the  horizontal  outpus  are  summed  and  then  a 
filter  is  applied,  has  been  treated  explicitly.  This  is  because  an  arbitrarily 
weighted  sum  of  2  horizontal  components  of  ground  motion  must  always  be 
equivalent  to  a  single  component  in  a  direction  defined  by  the  weights. 

Some  general  conclusions  may  be  stated  regarding  the 
usefulness  of  a  single  horizontal -component  instrument  for  eliminating 
noise  from  the  output  of  a  /ertical  instrument  at  the  same  location.  It 
would  appear  that  such  effects  as  system  noise,  gain  fluctuations  and  Love 
waves  need  not  present  serious  problems,  although  in  some  cases  they 
may  very  well  do  so.  The  most  important  consideration  is  usually  the 
properties  of  the  noise  to  be  eliminated.  If  significant  noise  rejection  is 
to  be  achieved,  then  some  rather  stringent  conditions  must  be  satisfied. 
Even  when  there  is  only  a  single-noise  mode  present,  it  is  essential  that  the 
azimuths  of  the  noise  sources  be  confined  within  a  range  somewhat  less 
than  180°  and  that  the  horizontal  be  oriented  near  the  center  of  th;s  range. 
This  requirement  follows  from  the  obvious  fact  that  a  filter  which  cancels 
out  noise  traveling  in  a  given  direction  must  also  amplify  noise  traveling  in 
the  opposite  direction.  In  practice,  such  extreme  directional  properties 
are  seldom  found. 
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If  more  than  one  noise  mode  is  important,  large  diff  ,rences 
in  the  horizontal -to-vertical  transfer  functions,  k(t),  cannot  be  permitted. 

In  Section  VII,  a  two-mode  situation  is  presented  in  which  no  noise  rejection 
is  possible  because  one  mode  is  characterized  by  a  prograde  particle  motion 
and  the  other  is  associated  with  retrograde  motion. 

With  a  3 -component  seismometer  at  point  location,  it  is 
theoretically  possible  to  determine  both  the  direction  of  propagation  and 
the  shape  of  the  particle  orbit  for  an  observed  Rayleigh  wave.  However, 
the  sense  of  propagation  and  the  sense  of  the  orbital  motion  cannot  be  found 
unless  a  spatial  separation  is  introduced  into  the  system.  For  this  reason 
it  is  concluded  that  ordinary  3 -component  systems  offer  little  promise  for 
most  applications  and  that  attention  must  be  directed  toward  arrays  of 
horizontal-  and  vertical -component  seismometers. 
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PART  II 


MULTICOMPONENT  ARRAYS 
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SECTION  IX 


FORMULAS  FOR  ARRAY  RESPONSE  IN  SINGLE -MODE 
NOISE  WITH  ARBITRARY  AZIMUTHAL  POWER  DISTRIBUTION 

A.  DESCRIPTION  OF  ARRAY 

Let  O  be  the  origin  of  an  X-Y  coordinate  system  in  the  plane 

of  the  earth's  surface.  We  consider  an  ai  ray  of  seismometers  consisting 

of  a  vertical  seismometer  v  located  at  O,  and  M  horizontal  seismometers 

h  ,  h_,  .  .  .  ,  h  ,  where  the  m^  horizontal  seismometer  h  is  located 
lc  M  m 

at  the  point  P  having  polar  coordinates  (r  ,  p  ),  and  h  is  oriented  in 
m  mm  m 

the  direction  tjt  (Figure  IX -1). 
m 


Y 


Figure  IX-1.  Generalized  Seismometer  Array 
Configuration 
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B.  ASSUMPTIONS  ABOUT  THE  NOISE  FIELD 


We  wish  to  examine,  for  this  array,  the  extent  to  which 
noise  on  the  vertical  can  b«e  predicted  from  noise  on  the  horizontals  in 
single -mode  trapped  noise.  We  make  the  same  assumptions  about  the 
noise  field  as  we  did  in  Section  III,  plus  two  additional  assumptions, 
stated  below.  To  simplify  the  calculations,  let  us  start  with  a  continuous 
azimuthal  power  distribution  function  P(8),  bearing  in  mind  that  this  is  the 
limiting  case  of  a  discrete  azimuthal  hower  distribution  function.  Recalling 
our  previous  notation,  we  have 

P(0)  »  azimuthal  power  distribution  function  (9.  1) 

1(f)  a  inline  horizontal  autopower  spectrum 

K(f)  =  inline  horizontal -vertical  transfer  function 

As  mentioned  above,  we  will  make  two  additional  assumptions 
about  the  noise  field;  these  assumptions  are  not  stringent,  and  they  are  made 
in  order  to  allow  us  to  derive  simple  expressions  for  the  crosspower  spectra 
between  spatially  separated  seismometers.  We  assume 

(1)  Noise  from  each  direction  propagates  in  plane  waves 

(2)  Let  h  and  h1  be  the  outputs  of  two  horizontal  seismometers 
located  at  points  A  and  B  respectively;  let  r  be  the  vector 
beginning  at  A  and  ending  at  B;  let  both  seismometers  be 
pointing  in  the  direction  r  .  Then  for  a  single  noise  source 
generating  waves  traveling  in  the  directior  r,  we  have 

h'  *  q®h,  the  Fourier  transform  of  q(t)  being  (9.2) 

Q(f)  *  exp  (2tt  i  r  k(f) ),  where  r  a  |r  |,  (9.3) 

k  a  wave  number 

[See  Laster  et  alj 
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Assumption  (1)  requires  that  the  distance  to  the  noise  sources 
not  be  small  relative  to  the  array  dimensions,  while  assumption  (2)  demands 
that  inline  horizontal  motions  at  two  separated  points  be  related  by  a  transfer 
function  Q(f)  which  is  pure  phase  shift  at  all  frequencies. 

C.  CALCULATION  OF  AUTOPCWER  AND  CROSSPOWER  SPECTRA 

Returning  to  our  array  of  seismometers  v,  h^ . h  , 

leu  us  derive  formulas  for 

i  ,  $  ,  etc. 

mn  vrn 


Let  us  first  compute  the  crosspower  spectrum  between  h^  and 

h  ,  horizontal  seismometers  located  at  P  =  (r  ,  p  )  and  P  =  (r  ,  p  ) 
n  m  mm  n  n  n 

respectively,  and  having  orientations  and  f  respectively. 

m  n 


Let  9  be  an  angle,  and  let  X  and  X  be  the  components  of 

m  n 

horizontal  motion  in  the  direction  0  at  P  and  P  respectively,  due  to  noise 

m  n  - - - - 

from  a  single  source  at  azimuth  9,  and  temporarily  let  h  and  h  denote 
- a -  m  n 

the  contribution  from  this  single  source  to  the  seismometer  outputs. 


Then 

h  -  -cos  (9  -  \|f  )  X 
m  mm 

h  =  -cos  (9  -  'll  )  X 
n  n  n 

and,  by  assumptions  (1)  and  (2), 

X  =  q  ©  X 
m  n 

where  the  Fourier  transform  of  q(t)  is 

Q(f)  =  exp  [2™  i  k(f)  (r  cos(0  -  P  )  -  r  cos(6  -  P  )  )| 

1  m  m  n  n  1 


(9.4) 


(9.  5) 


(9.6) 
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Therefore,  for  a  single  wave  from  direction  0,  we  have,  by 
(?.  4)  -  (9.  6)  that 


$  =  cos  (0  -  t  )  cos(0  -  l|f  )  . 

mn  m  n  a.  k 

m  n 


(9.7) 


=  cos  (0  -  \|f  )  cos  (0  -  'll  )  Q  $. 

m  n  A.  a. 


n  n 


cos(0  -  til  )  cos(0  -  >li  )  exp  Ett  i  k  (r  cos(0  -  p  ) 
m  n  m  m 


-rncos(0-  Pn))  |  ix 


n  n 


[Crosspower  spectrum  for 

single  source 


To  obtain  the  crosspower  spectrum  for  the  sum  of  nc*  •  f  .  ".i 

all  sources,  we  set  ^  =  P(0)#$  in  (9.7)  and  integrav.  over  all  values  of 
0  from  -n  to  TT.  n  n 

Therefore  the  crosspower  spectrum  between  h  ar  h  in  a 

m  n 

single -mode  noise  field  with  azimuthal  power  distribution  function  P(0)  is 

f  n 

$  =  $  I  cos(0  -  'll  )  cos(0  -  'll  )  exp  Iett  i  k  (r  cos(0  -p  )  (9.  8) 

mn  -,_tt  m  n  '  m  m 

-r  cos(8  -  D  ))|  P(0)  d0 
n  n 

In  case  m  =  n,  (9.  8)  reduces  to 


f  n 

i  =  $  cos2  (0  -  ’|i  )  P(0>  d  0 
mm  -tt  m 


(9.9) 


Cf.  Eq.  3.  10 
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A  similar  computation  yields 


f  n 

$vr="K$i  cos  (S  -  )  exp  |2TT  i  k  r  cos^-P  )|P(Q)dc  (9.10) 

vm  •Ltt  m  1  m  rn  I  ' 

where  K (f)  is  the  horizontal-vertical  transfer  function.  (See  Eq.  9.  1). 


.  As  before  (Eq.  3.  12), 


$  =  Ik  |  * 

vv 


(9.11) 


D.  PREDICTION  FILTERS 

The  optimum  prediction  filter  transfer  functions  G  , 

m 

m  =  1,  2,  .  .  .  ,  M,  are  determined  by  the  system  of  equations  (Burg,  1964) 


'•11 

$  $ 

12  *  ‘  ‘ 

V 

V 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

$ 

Ml  ' 

.  $ 

MM 

gm 

■»  <■ 

$ 

vM 

where  the  $  ,  $  are  given  in  Equations  (9.  8)  -  (9.  10). 

mn  vm  0  ' 
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SECTION  X 


ISOTROPIC  SINGLE -MODE  NOISE 


In  this  section  we  shall  consider  the  theoretical  capability  of 
particular  arrays  of  spatially  separated  horizontal  seismometers  for  pre¬ 
dicting  a  vertical  component  of  single -mode  isotropic  noise.  Array  geo¬ 
metries  will  be  specified  by  use  of  the  parameters  r  ,  p  ,  iji  ,  and 

m  m  m 

use  will  be  made  of  formulas  (9.  8)  -  (9.  12). 

A.  INTEGRAL  REPRESENTATIONS  OF  BESSEL  FUNCTIONS 


In  order  to  simplify  the  expressions  which  will  be  obtained 
in  this  section,  the  following  integrals  are  useful.  Watson  (p.  41)  gives  the 
following  form  for  Poisson's  integral  representation  of  the  Bessel  function 
of  first  kind  and  of  order  v. 

J  (z)  -  — | - - —  f  cos(z  cos  6)  sin^  ;  6  d8  (10.  1) 

v  r(v+y)  r(j)  0 

Noting  t.iat 

r(y)  =  ^  andr(3/2)  =  A/ 2, 

J  (z)  =  —  f  cos(z  cos  8)  d8  =  r— —  f  cos(z  cos  0)  d9  ^ 
0  tt  ^  2  n  J  ^  ' 


2  CU 

"  Jo 


.  TT 

2  If  2 

cos  (z  cos  0)sin  0d0=£^J  cos(z  cos  0)  sin  0  d  0 


(10.3) 


-TT 


Webster  (p.  322)  gives,  the  recurrence  relation 

J,U>  , 

—  =  J  Jj0(.)t  J2(z)j 
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Hence, 


B.  2 -CHANNEL  SYSTEM:  SEPARATED  VERTICAL  AND  HORIZONTAL 
Recall  from  Section  III  that  no  prediction  of  a  -'ertical 
component  from  a  horizontal  component  located  at  the  same  point  is  possible 
in  isotropic  noise.  This  fact  is  one  of  the  most  severe  limitations  on  the 
potential  use  of  two -component  seismometers  at  point  location,  since  the 
noise  at  many  recording  stations  appears  to  be  nearly  isotropic. 

However,  a  portion  of  noise  on  a  vertical  seismometer  can 
be  predicted  from  a  spatially  separated  horizontal  seismometer,  even  in 
isotropic  noise;  it  is  this  fact  that  motivates  the  study  of  arrays  of  the 
general  type  discussed  in  the  last  section. 

Let  us  begin  by  examining  the  prediction  capabiKtv  of  the 
following  system  in  single -mode  isotropic  loise: 
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SYSTEM  0;  Vertical  seismometer  v  located  at  the  origin  O,  and 
a  single  horizontal  seismometer  h  at  distance  r  from  the  origin,  with 
azimuth  p  =  0,  and  orientation  ilr.  (Figure  X-l) 

Since  the  noise  is  isotropic,  we  substitute  P(9)  =  —  in 

2tt 

formula  (9.  10)  to  obtain  $  ,  .  We  get 

vh  ° 


-K  $ 

2  TT 


cos  (9  - 


i2n  rk  cos  9  , « 
e  do 


(10.7) 


-iK  i 
2H  ' 


f  n 

cos  9  cos  ilr  sin  (2tt  rk  cos  9)  dG 


=  -iK$  •  cos  \|f  *  Jj  (2n  rk) 

Mow,  by  (9.  9)  and  (9.  10)  we  have 

tt 

#hh  =  iJ  c°*2 <e •  *• de  *  j  *  <10-8> 

-TT 


Therefore  the  coherence  between  v  and  h  for  system  0  in  a 
single -mode  isotropic  noise  is 


Figure  X-l.  Geometry  of  System  0 
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Notice  in  (10.9)  the  dependence  of  coherence  on  'f  ,  the 

TT 

horizontal  seismometer  orientation  angle.  For  iff  =  ±—  ,  we  have  C  =  0; 

Ml 

and  for  f  =  0,  11  we  have  the  maximum  coherence, 


C  = 


(10.  10) 


Let  us  give  the  name  System  1  to  the  special  instance  of 
System  0  when  i|f  =  0  (Figure  X-2).  In  Table  X-l,  the  prediction  error 
1-C  for  System  1,  where  C  is  calculated  by  (10.8),  is  given  as  a  function 
of  2tt  rk.  With  this  two-channel  system,  we  get  a  prediction  error  as  low 
as  0.  324,  for  2tt  rk  =1.8. 


C.  MULTIPLE -HORIZONTAL  ARRAYS 


Since  par1-  of  the  vertical  component  of  single -mode  isotropic 
noise  can  be  predicted  by  a  single  spatially  separated  horizontal  seismometer, 
one  should  expect  that  ver;f  good  performance  might  be  achieved  by  imple¬ 
menting  multichannel  prediction  using  a  large  enough  number  of  additional 
horizontals,  arranged  in  a  suitable  geometry.  To  nvestigate  this  possibility, 
let  us  consider  the  arrays  shown  in  Figure  X-2.  The  array  parameters  are: 


SYSTEM  1. 

V 

rl 

= 

r, 

pi 

=  t.  =  0 

SYSTEM  2. 

h  : 

r 

— 

r, 

p 

=  4  =  Ini 

m 

m 

m 

m 

SYSTEM  3. 

h  : 

r 

— 

r, 

p 

•  ♦  =  (m:' 

m 

m 

m 

m  3 

SYSTEM  4. 

h  : 

r 

r, 

P 

=  t  =  mTT , 

m 

m 

m 

m 

SYSTEM  5. 

h  : 

r 

= 

r, 

P 

M  =2^ 

m 

m 

m 

m  2 

SYSTEM  6. 

h  : 

r 

= 

r, 

P 

=  A  =  2HL 

m 

m 

m 

m  3 

(1-H) 


,  m  =  1,  2,  3 


,  m  =  1 ,  2,  3,  4 


SYSTEM  7.  h  :  r  =  r,  p  =  I } 
-  mm  mm 


2mTT 


,  m  =  1, 2,  3 


In  all  of  these  systems,  the  horizontals  are  place  at  the  same 
distance,  r,  from  the  vertical  seismometer  location,  and  the  horizontals  are 
oriented  so  that  they  point  away  from  the  vertical,  i.  e.  ,  f  =  p. 


X-4 


SYSTEM  1 


SYSTEM  2 


SYSTEM  3  SYSTEM  4 


\  / 


Figure  X-2.  Array  Geometries 


Table  X-l 


PREDICTION  ERROR  FOR  SYSTEM  1 


Dimensionless  Wavenumber,  Prediction  Error 

2rr  r  -k  (1-C) 

(radians) 


0 

0.  2 
0.4 
0.  6 
0.  8 
1.0 
1.  2 

1.4 
1.6 
1.8 
2.0 
2.  2 

2.  4 
2.  6 
2.  8 
3.0 

3.  2 

3.4 
3.  6 
3.8 
4.0 


1.0 
0.  980 
0.  923 
0.836 
0.728 
0.613 
0.  503 
0.413 
0.  350 
0.  324 
0.  335 
0.  382 
0.459 
0.  557 
0.  664 
0.770 
0.  863 
0.936 
0.  982 
1.000 
0.991 
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The  decision  to  restrict  our  attention  to  such  outward¬ 
pointing  arrays  is  motivated  by  the  following  fact:  let  us  be  given  a  set  of 

points  P  ,  P  ,  .  .  .  ,  P  in  the  X-Y  plane.  Then  among  all  arrays 
1  2  M 

V,  h  ,  h_,  .  .  .  ,  h  such  that  the  vertical  v  is  located  atO  and  the 
1  M 

horizontal?,  h  are  located  at  P  ,  m  =  1 ,  2,  .  .  .  ,  M,  with  h  having 
mm  m 

arbitrary  orientation  angle  ,  the  array  which  gives  the  best  multi  - 

m 

channel  prediction  performance  in  single-mode  isotropic  noise  is  the 

outward-pointing  array,  that  is,  the  array  for  which  lit  =  p  for  all  m. 

'  mm 

This  follows  from  the  fact  that  maximum  coherence  between  separated 
vertical  and  horizontal  in  single -mode  isotropic  noise  occurs  when  the 
horizontal  is  outward -pointing  (Formula  10.9). 

Having  restricted  ourselves  to  outward -pointing  arrays, 
we  further  restrict  this  preliminary  investigation  to  arrays  in  which 
all  of  the  horizontals  are  at  equal  distance  from  the  vertical  because,  for 
an  array  of  outward -pointing  horizontals  h  at  equal  distance  from  the 
vertical  v,  the  theoretical  optimum  prediction  filters  to  predict 
v  from  the  h^'s  in  isotropic  noise  are  quite  easy  to  compute.  In  fact, 
for  each  of  the  Systems  1,  2,  4,  5,  6,  and  7  it  is  obvious  that  all  of  the 
prediction  filters  are  equal,  because  of  symmetry  both  in  the  noise  field  and 
in  the  array  geometries.  Hence,  it  is  necessary  to  design  but  one  prediction 
filter  for  each  of  these  arrays,  and  in  each  case  apply  this  filter  to  the 
summed  outputs  of  all  the  horizontals.  Therefore,  the  minimal  prediction 
error  for  these  systems  is  simply 

Prediction  error  =  1 -C  (10.12) 

vs 

where  s  =  summed  output  of  horizontals,  and  C  =  coherence  between  v  and  s. 

vs 
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One  can  obtain  expressions  for  C  in  terms  of  Eessel 

vs 

functions  for  each  of  the  systems  1-7.  In  the  case  of  system  3,  1  -Cvg 
may  not  be  the  minimal  prediction  error,  since  it  is  possible  that  the 
optimum  prediction  filters  for  h  ,  h  ,  h  are  not  identical. 

U  1  w 

In  the  case  of  radially -oriented  horizontal -component 
instruments  on  a  circle  of  radius  r  in  isotropic  noise  (with  total  power  $), 
the  expressions  (9.  8)  and  (9.  10)  may  be  simplified,  as  follows: 

In  (9.8)  let 

<j>+a  =  0  -  ib  =  0  -  p 

mn  n  n 

<i>  -  a  =  0  -  =  0  -  p 

mn  m  m 

u  =  2nk  r 

P(0)  =  $  /2tt 


Thus  the  crosspcwer  between  two  horizontals  is 

$  r  n 


,  *  f  f  2  .  2  .  .  2  n  2  1 

*  =  I  <  cos  a  sin  <i>  -  sin  a  cos  <b> 

mn  £iT  j  ^  ^  mn  mn  ) 


cos(2u  sin  a  cos<(>)  dcj> 
mn 


=  $  |cos(2a  )  J  (2u  sin  a  )  +  J  (2u  sin  a  )}/2 
t  mn  U  mn  c  mn  J  j 

where  2  a  =  p  -  p  is  the  angular  separation  between  instruments  m 
mn  n  m 

and  n. 

The  crosspower  between  the  central  vertical  and  a 
horizontal  is  given  by  substituting  in  (9.  10): 


d,  =  0  -  lb  =  0  -  p 

m  m 

u  =  2tt  kr, 


(10.  1 
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giving 


$ 

vm 


-iK  $ 
2tt 


cos  <|>  sin(u  cos4>) 


=  -i  K  $  •  J  (u) 


The  autopower  for  a  horizontal  is  $  =  $/2.  This 

mm 

mav  be  verified  by  setting  CX  =  0  in  (10.  13).  The  autopower  for  the 

mn  r 

vertical  is  =  |Kp  $  as  given  by  (10.8). 


The  coherence  between  the  output  of  the  vertical  and  the 
summed  outputs  of  the  horizontals  is 

$  .  $ 

£  _  vs _ sv 

vs  $  $ 

vv  ss 


where 


and 


$ 

vs 


M 

=  V  $ 

vm 


=  -  i  M  K  $  J^u) 
M  M 

$  =  y.  yi  $ 

s  s  ^  mn 

m=l  n=  1 


=  $ 


m=  1  n  =  l  { 


cos  (2a  )  Jn(2u  sin  a  ) 

mn  u  mn 


J_(2u  sin  a  )  i  /2 
2  mn  )  f 


for  an  array  containing  M  horizontals. 


(10.  14) 


(10.  15) 


(10.  16) 


(10.  17) 
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Thus,  ,  2  (10.18) 

r  .  2M  |Ji(u>| 

vs  MM 

EV"  icos  (2a  )  J  (2u  sin  a  )  +  J  (2u  sina  )i 

i—i  (  mn  0  mn  2  mn  ) 

m~  i  n=i 

Denoting  the  coherence  for  system  1  by  C*  ,  -he  coherence  for  system 
2  vs  7 

2  by  C  ,  etc.  ,  we  have: 
vs 

cts  *  iji(u,r/({)  <ioi9> 

Cvs  =  4|J1(u)|2/{‘  +  J2 

cvs  *  9K<u>|2/{1  +  J° <u)  J0  ('/3_u)  +  2J2'u»  +  J2(^u)} 

Cvs  '  *W7  i1  -  V2U>  +  J2(2u»} 

CVS  =  16jj,(u)|2/{2+4J2(#-u)  +  2J2  (2u)  -  2J0  (2u)} 

Cvs  =  36|Jl(u)|  /{3  +  3  J0  <“>  +  6  J2  <u>  -  3J0  (^u) 

+  6  J2  (y/lu)  +  3  J2  (2u)  -  3J0  (2u)| 

Cv.  -  +  3  J2^“)  -  1  J0  (/*"»)} 

The  prediction  errors  (1-C^)  for  the  seven  systems  are 

illustrated  in  Figures  X-3  through  X-9  by  the  curves  labeled  "P.  E.  " 

In  addition  to  the  prediction  error,  there  are  certain  functions  which 

determine  the  usefulness  of  each  array.  It  is  not  sufficient  that  the 

prediction  error  be  low,  if  the  power  response  $  is  also  low,  since 

s  s 

the  optimum  prediction  filter  response  is  inversely  proportional  to 
and  it  is  desirable  to  avoid  undue  amplification  of  uncorrelated  noise 
appearing  on  the  horizonal  outputs. 
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Response  F ^motions  for  System 


Response  Functions  for  System  3 


Response  Functions  for  System  4 
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Figure  X-7.  Response  Functions  for  System 


SYSTEM 


Response  Functions  for  System  7 


A  measure  of  the  amplification  of  Rayleigh  noise  relative  to 

uncorrelated  noise  is  shown  by  the  curves  marked  "R"  in  Figures  X-3 

through  X-9.  The  function  plotted  is  the  ratio  of  Rayleigh  noise  power  in 

2 

the  average  of  the  M  horizontal  outputs  )  to  the  Rayleigh  noise 

s  s 

power  in  the  output  of  a  single  horizontal  ($/2).  The  corresponding  ratio 
for  uncorrelated  noise  is  1/M. 


D.  LOVE -WAVE  RESPONSES 


Another  important  consideration  is  the  response  of  the 
system  to  Love  waves.  Since  the  vertical  instrument  is  insensitive  to 
Love  waves,  any  Love  energy  which  appears  in  the  average  of  the  hori¬ 
zontals  will  lead  to  erroneous  prediction  of  noise  on  the  vertical.  It  is 
therefore  desirable  that  the  system  be  as  insensitive  as  possible  to 
Love  waves.  In  Figures  X-3  through  X-9,  the  curves  labeled  "L" 

show  the  ratio  of  Love  noise  power  in  the  a/erage  oi  the  M  horizontal 

2 

outputs  (  A  /M  )  to  the  Love  noise  power  in  the  output  of  a  single 
s  s 

horizontal. 


The  Love -wave  responses  are  computed  as  follows:  For  our 
purposes,  we  may  consider  a  Love -wave  to  be  identical  to  a  Ravleigh  wave 
except  that  the  motion  is  perpendicular  to  (instead  of  inline  with)  the  direction 
of  propagation  and  there  is  no  vertical  component  of  motion.  Assume  a 
uniform  power  distribution  A/2tt.  Since  there  is  no  vertical  component, 


A 


vm 


vv 


=  0. 


Since  the  horizontal  motion  is  now  transverse,  the  factors 

cos(9  -  !jf  )  and  cos  (9  -  i|f  )  in  equation  (9.  8)  must  be  replaced  by 
m  n 

sin(0  -  t  )  and  sin(0  -  ijf  ).  By  means  of  the  same  substitutions  as 
m  n 

were  used  to  obtain  Equation  (10.  13),  we  derive  the  analogous  expression 
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A  f  <  2  2  .2  2  i 

A^„  =  T7T  ccs  a  cos  o  -  sin  a  sin  cos(2u  sin  a.  cosd>)  d<b 

mn  2n  J  n  (  mn  T  mn  Tf  '  mn  Y 

=  A  |cos  (2a  )  J  (2u  sina  )  -  J  (2u  sin  a  )l/2 

I  mn  0  mn  2  mn  1/ 

The  power  spectrum  for  the  sum  of  M  horizontals  is 

A  =  A  2_.  <  cos  (2a  )  J  (2u  sin  a  )  -  J  (2u  sin  a  )  i  /  2 

ss  *-*  (  mn  0  mn  2  mn  f  / 

m=  1  n=l 

Rayleigh -wave  response  functions  §  /$  for  the  7  systems 

are  given  by  the  denominators  in  Equations  (10.  19).  Corresponding  Love- 
wave  response  functions  A^/A  may  be  obtained  by  changing  the  signs 
of  all  terms  in  the  corresponding  Rayleigh  functions  which  contain 
Bessel  functions  of  order  2. 
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SECTION  XI 


CONCLUSIONS 

Theoretical  results  obtained  in  this  study  indicate  that  arrays 
of  horizontal-component  seismometers  should  prove  to  be  useful  tools  for 
the  removal  of  trapped-rnode  noise  from  the  outputs  of  vertical-component 
instruments.  A  necessary  condition  is  that  the  horizontals  must  be  spatially 
separated  from  the  vertical.  In  the  case  of  a  single -noise  mode,  arrays 
such  as  those  studied  in  Section  X  should  perform  best  when  the  array 
diam  r  is  approximately  one -half  wavelength.  Signal  enhancement 
systems  employing  only  vertical  instruments  require,  in  general,  array 
diameters  of  at  least  one  wavelength.  Thus,  multicomponent  arrays  should 
offer  meaningful  advantages  in  terms  of  land  and  telemetry  requirements. 

It  is  unlikely  that  difficulties  presented  by  system  noise, 
uncorrelated  seismic  noise  and  Love  waves  should  be  any  more  serious 
than  they  are  in  the  case  of  vertical -component  arrays.  The  effect  of 
additional  noise  modes  on  multicomponent  array  performance  has  not  been 
studied  yet.  However,  it  is  reasonable  to  assume  that  a  multiplicity  of 
modes  may  be  dealt  with  by  the  application  of  multichannel  filter  techniques 
to  the  outputs  of  rings  of  horizontals.  It  has  been  shown  that  an  array 
consisting  of  a  number  of  rings  of  verticals  can  be  useful  in  the  presence 
of  a  similar  number  of  noise  modes,  and  there  is  no  known  reason  for 
assuming  that  this  usefulness  might  be  a  property  peculiar  to  vertical - 
component  arrays. 
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Tne  results  derived  theoretically  for  multicomponent 
seismomter  arrays  are  quite  conducive  to  optimism.  Referring  to 
Figures  X-7  and  X-8,  it  can  be  seen  that,  if  the  array  dimensions  are 
suitable  (r/X~0.25),  the  theoretical  prediction  error  is  less  than  0.01, 
the  Love  wave  response  function  is  less  than  0.  01  and  the  Rayleigh  wave 
response  function  is  approximately  0.7.  Thus,  it  is  theoretically  possible 
to  predict  (and  hence  remove)  more  than  99  percent  of  the  Rayleigh  noise 
on  the  vertical  component  in  the  frequency  range  in  which  the  wavelength 
is  appropriate.  The  low  value  of  the  Love  wave  response  function  implies 
that  very  little  extraneous  noise  power  should  be  introduced  into  the 
system  output  as  a  result  of  Love  waves  appealing  in  the  outputs  of  the 
horizontal  instruments.  The  high  value  of  the  Rayleigh  wave  response 
function  implies  that  the  filter  responses  need  not  be  unduly  large  and 
that  uncorrelated  noise  from  the  horizontals  will  not  be  amplified  to  a 
serious  degree.  The  finding  that  all  three  response  functions  take  on 
desirable  values  in  the  same  range  of  wavelength  is  both  unexpected 
and  fortuitous. 

From  the  results  of  the  theoretical  investigations  reported 
in  this  report,  it  is  concluded  that  multicomponent  seismometer  arrays 
offer  a  great  deal  of  promise  for  signal  enhancement  applications.  It 
is  recommended  that  experime’- ' *•  ]  investigations  of  such  systems  be 
undertaken  as  soon  as  possible. 
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